Coherent Oscillations and Giant Edge Magnetoresistance in Singly Connected 

Topological Insulators 



o 
o 

U 

Or 
< 



i 

C 

o 
o 



> 

oo 

o 
o 

OS 

o 



X 



Rui-Lin Chu, 1 Jian Li, 2 J. K. Jain, 3 and Shun-Qing Shen 1 

'Department of Physics and Center for Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
2 Department of Theoretical Physics, University of Geneva, CH-1211 Geneva 4, Switzerland 
'Department of Physics, 104 Davey Lab, Pennsylvania State University, University Park, PA 16802, USA 

(Dated: April 6, 2009) 

The topological insulators have a gap in the bulk but extended states at the edge that can carry 
current. We study a geometry in which such edge states will manifest themselves through periodic 
oscillations in the magnetoconductance of a singly connected sample coupled to leads through narrow 
point contacts. The oscillations occur due to quantum interference of helical edge states of electrons 
traveling along the circumference of the sample, and have a period of AB — h/eA a g, where A e g is 
the effective area enclosed by the edge states of the sample. Our calculation indicates the possibility 
of a large change in the magnetoresistance at small B, termed giant edge magnetoresistance, which 
can have potential for application. The zero field conductance also exhibits oscillations as a function 
of the Fermi energy due to interference between edge states. The amplitude of the oscillations is 
governed by, and therefore can be used to investigate, the transverse width of the edge channels. 

PACS numbers: 73.23.-b, 73.43. Qt, 72.25.-b 



Topological insulators differ from ordinary insulators 
in having a pair of extended helical edge states, which 
results in quantum spin Hall (QSH) effect 0, 0, 0|. 
Several candidates for topological insulators have been 
proposed, and a non-zero conductance has been mea- 
sured experimentally in the "inverted" -band semicon- 
ductor HgTe/CdTe quantum well in a band insulat- 
ing region. *3] More experiments are beginning to ex- 
plore the edge and surface state properties in topological 
insulators. [J [f| H, 0, HI A direct observation of extended 
edge states would be important to establish the physics of 
topological insulators, and one may ask if they can also 
exhibit other interesting phenomena. We demonstrate 
that an effective one-dimensional ring is formed between 
two consecutive scatterers, which leads to Aharonov- 
Bohm (AB) oscillations in conductance. An observation 
of such oscillations in a singly connected geometry will 
constitute a direct observation of edge transport. Such 
oscillations are analogous to similar oscillations in singly 
connected quantum Hall systems 

HEnm. but with an 

important difference: in the present case, the oscillations 
occur at very small magnetic fields. Our results also indi- 
cate the possibility of a "giant edge magnetoresistance" 
(GEMR), which is insensitive to the geometry of the de- 
vice and may have potential for practical application. 

We study here the device shown in FIGQ] which con- 
sists of a two-dimensional strip of a topological insula- 
tor on which two quantum point contacts (QPCs) have 
been patterned thorough gates (shaded regions in FIG. 
[T]). The QPCs define a saddle shaped confining potential, 
whose height can be controlled by a gate voltage. An ef- 
fective disk of area A — ttR 2 (R is the radius of the disk) 
is formed in the center. An AB effect in the device can 
be expected intuitively because a topological insulator 



FIG. 1: Schematic of the geometry studied in this manuscript, 
which consists of a disk connected to two reservoirs through 
two quantum point contacts. Red and blue lines indicate the 
chiral edge channels of spin up and down electrons, with ar- 
rows indicating the direction of their motion. In our calcula- 
tions, we take the following parameters: width W — 500nm, 
length L — 1030nm; the slit width d of the quantum point 
contacts is taken as a variable. 



possesses a pair of independent gapless edge states of dif- 
ferent spins moving in opposite directions, each forming 
an ideal one-dimensional loop around the disk. The two 
edge states are independent because no backscattering is 
allowed at a given sample edge even in the presence of 
weak time-reversal invariant disorder. We note here that 
spin is not a good quantum number in topological insu- 
lators because of spin orbit coupling. In the absence of 
a magnetic field, the actual edge states are eigenstates of 
the time reversal operator; their characterization as spin 
up and down is not precisely correct, and the word "spin" 
below is to be viewed more generally as the quantum 
number denoting the two states of a Kramers doublet. 
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The effect of magnetic field, which breaks time reversal 
symmetry, is complicated. In our considerations below 
we assume that the magnetic field is sufficiently weak 
that the Zeeman term can be neglected; in that limit, 
the field does not open any gap in the spectra for the 
edge states, as discussed in detail in Ref. [l 21 ] . 

In the remainder of the article we study this geometry 
quantitatively. We will assume that the weak magnetic 
field B is normal to the plane. Following Ref. [l0( we 
consider a spin-up (or spin-down) electron travelling from 
the left hand side (LHS) in FIG. 1. At the LHS junc- 
tion it splits into two partial waves: one is transmitted 
through the QPC into the disk with amplitude t, and the 
other is transmitted across the QPC with an amplitude 
r causing a backscattering. We denote the wave func- 
tion amplitudes in the upper edge and lower edge, right 
after the LHS junction, by u\ and d\ 1 respectively. The 
corresponding amplitudes in the vicinity of RHS junc- 
tion are 112 = u\ exp[i#/2] and d,2 = d\ exp[— i6/2], where 
6 = 2n(j)/ 4> a + 86, 4> a = h/e a fundamental unit of mag- 
netic flux, (j> = ttR 2 s B is the magnetic flux threading 
the effective one-dimensional loop with an effective ra- 
dius Reg, and 86 is the phase acquired by the wave func- 
tion traveling along the loop 86 — 2 / KkR e g- A partial 
wave goes through the RHS slit with an amplitude t' 
and across the slit with an amplitude r'. To simplify the 
problem, we assume identical reflection and transmission 
amplitudes for the two slits, t = t' and r = r'. Using the 
theory of multi-scattering processes [l3| it then follows 
that the total transmission for spin-up electron through 
the slit A and B is given by 



T T (B) = 



To 2 +4(1 -To) cos 2 



(1) 



Here To = 1 1 2 | is the transmission coefficient of an elec- 
tron through a slit, and #0 = arg(rr'). Resonant tun- 
neling occurs for cos 6+ ®° = 0, i.e., T^(B) = 1 for any 
T . The transmission coefficient for a spin-down electron 
T^(B), which is the time reversal counterpart of spin- 
up electron at —B field, is given by T±(B) = TU-B). 
According to the Landauer-Buttiker formula, fl4 . 15 1 the 
total conductance is 



G(B) = — [T^(B)+T l (B)}. 



h 



(2) 



The AB oscillations in the conductance G as a function 
of the magnetic flux (f> through the disk are therefore 
expected to be symmetric with respect to the direction 
of the magnetic field. 

To make further progress, we undertake a numerical 
calculation for the specific case of an "inverted" -band 
HgTe/CdTe heterojunction. We consider with the effec- 
tive Hamiltonian for HgTe/CdTe quantum wellQ 




/ h(k) 
n \ h*(-k) 



(3) 



E(meV) 



FIG. 2: (a) Magnetoconductance G for several energies in 
the band gap. (b) The locations of the energies for the curves 
in panel (a), (c) Conductance G T as a function of the Fermi 
energy at zero magnetic field (B = 0). All calculations assume 
point contact constriction width of d = 40nm. 



where, in the spin-up sector, h(k) = —Dk 2 + A(k x a x + 
kydy) + (M — Bk )<j z , and k = (k x ,k y ) is the wave vector 
in two dimension, k is the Pauli matrix vector. A, B, D 
and M are sample specific parameters, and are functions 
of the thickness of the quantum well; here we take A — 
364.5meV nm 2 , B = — 686meVnm 2 , D — — 512meVnm 2 , 
AI = — lOmeV. In the spin-down sector h*(—k) is the 
time reversal counterpart of h(k). The existence of edge 
states in this model has been discussed previously. [13, Ha] 
To determine the conductance through the disk sub- 
jected to a perpendicular and weak field B, we use 
the Keldysh Green function technique to calculate the 
transmission coefficients T numerically in the Landauer- 
Biittiker formalism [13|. For this purpose, we use a tight- 
binding model that reproduces the Hamiltonian in Eq. 
(3) in the continuum limit, and include the magnetic field 
through the substitution k — > k — f A (A is the vector 
potential). In the present work, more 400 x 200 lattice 
sites and the lattice space a — 2.5nm were used in the 
calculation, and the lattice size effects are vanishingly 
small, so all results reported here reflect the continuum 
limit. 

FIG. 2 shows the conductance for several Fermi ener- 
gies, whose locations in the band structure are illustrated 
in FIG. 2(b). The conductance exhibits periodical oscil- 
lations in B, with the period determined by the magnetic 
flux. For most parameters, the maximal value of the con- 
ductance does not reach 2e 2 /h because T^{E) and T' L (B) 
in general do not satisfy the resonance condition simulta- 
neously. A new feature is the rapid oscillations in the con- 
ductance at B = as a function the Fermi energy Ef, as 
shown in FIG. 2(c). These oscillations are caused by the 
phase shift 2irkR e s wher e the dispers ion of the edge state 
E(k) ~ -MB IB ± Ay/1 -D 2 /B 2 k w +7.46 ± 242.5A; 
meV (k in unit of nm -1 ) is linear in k. The increas- 
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FIG. 3: The magneto-conductance of a single spin channel 
for several widths d of the quantum point contacts, obtained 
numerically for the sample of Fig. 1. The black curves are 
fits using Eq. |TJ. Smaller d produces smaller t (CHECK) 
and hence larger oscillation amplitude. The insert shows the 
relation between To and d. The period is independent of d. 
All curves are evaluated for Ef — 0.5meV. 



FIG. 4: Panels (a) and (b) show the local current density 
distributions on the sample of Fig.l for magnetic fields B 
at a peak and a valley of the magnetoresistance, which cor- 
respond to maximum transmission and maximum reflection, 
respectively. The current reflection is totally suppressed in 
(a). The panels (c) and (d) show the corresponding local 
density of states distributions. The calculation assumes the 
parameters d = 30nm and E — 0.5meV. 



ing oscillation amplitude reflects that the spatial dis- 
tribution of edge state is dependent on the Fermi en- 
ergy. Accordng to Zhou et al. 12] the spatial distri- 
bution of the edge state near the boundary has the 

form (j>{y) = {e~ v /^ - e^/^)/^/% + % - |^ where 

d « 317/(240 + E) and £ 2 w 317/(13.4 - E) nm (E 
in units of meV) near the crossing point k = for the 
parameters adopted in this paper. Because we typically 
have £ x << the spatial distribution of the edge state 
is determined dominantly by £ 2 . Thus the edge state be- 
comes wider with increasing energy, ultimately evolving 
into a bulk state when E > 13.4meV. Correspondingly, 
the transmission coefficients decreases and the oscillation 
amplitude increases according to Eq.(l). 

FIG. 3 shows the dependence of oscillations on the slit 
width d. The solid lines are the best fits from Eq.(l) with 
three adjustable parameters To, i? ff and 59. The effec- 
tive radius i? ff = -\/20 o /i?o can be deduced from the 
period Bo of the oscillation; its ratio to the radius of the 
disk is R e s/R ~ 0.95, which reflects the finite width of 
an edge state near the boundary of the disk. The inset in 
FIG. 3 shows how T depends on the slit width d. As ex- 
pected, To vanishes for d/^ 2 << 1 an d approaches unity 
for d/£ 2 >> 1- Conductance oscillations are suppressed 
in both limits; in particular, for d/£ 2 >> 1 it has the 
universal value 2e 2 /h [3, 0] and displays no magnetore- 
sistance. The AB oscillations occur only when the slit 
width is comparable to the spatial distribution of edge 
state. The numerical results are in good agreement with 
Eq.(l) in the first one or two periods, but exhibit an 
increasing deviation for a larger field B, which may orig- 
inate from the neglect, in our fits, of the dependence of 
the energy dispersion of the edge states on magnetic field 



(we use the B — dispersion when fitting the numerical 
data). 

The quasi-one dimensional quantum interference be- 
havior can be further illustrated by plotting the spatial 
profile of local current and the local density of states 
in a nonequilibrium situation [17[, where Fermi surface is 
slightly higher in the left lead so that the electrons flow 
from left to right. Current flow patterns for the up spin 
sector are shown in Figs. 4(a) and 4(b). The current 
flows only along the boundary of the device with defi- 
nite chirality (see FIGQ]). At the resonance transmission 
To = 1, the backward current reflection at QPC A is 
totally suppressed, which corresponds to a peak in the 
magnetoresistance oscillations. The maximum reflection 
is observed at the valley of the oscillations. When a reso- 
nance transmission occurs, the maximum of local density 
of states is found in the central disk edge while, at mini- 
mum transmission, the maximum local density of states 
is observed in the terminals. 

The feasibility of the observations of the predicted 
magnetoresistance oscillations depends in part on how 
sensitive their amplitude is to the ever present disor- 
der. Given that time reversal invariant disorder does not 
cause any backscattering between two spin channels at 
the same edge[H [2(1 21 1, one might expect that it also 
does not affect the coherent magnetoresistance oscilla- 
tions. We study the effect of disorder by introducing 
disorder as a random on-site energy with a uniform dis- 
tribution within [— w/2, w/2]. The results, displayed in 
FIG. 5(a), demonstrate that disorder diminishes the am- 
plitude of oscillations while enhancing the average con- 
ductance. The reduction of the coherent oscillations orig- 



4 




FIG. 5: (a). Conductance oscillations in the spin up sector 
as a function of disorder. The oscillations gradually flatten 
but electron tunnelling through the slit is significantly en- 
hanced with increasing disorder strength, (b) The disorder 
dependence of transmission coefficient through a single QPC 
To for three different spatial distributions of disorders. Red 
line - disorder is distributed over the entire device. Blue line 
- disorder is distributed only in the vicinity of the QPC slit. 
Pink line -no disoder in the vicinity of the QPC slit. The 
inset in (b) shows the variance of the period of magnetic field 
as a function of the disorder strength for two different spa- 
tial distribution of disorder. Averaging is performed over 200 
disorder samples; the parameters are chosen as d — 40nm, 
E = 5.5meV. 



inates from an impurity configuration dependent phase 
shift 89 in Eq. (1); the disorder averaging over the phase 
shift smears the quantum interference, completely sup- 
pressing it when the variance of 89 is comparable with n. 
It is noteworthy that the AB oscillations survive up to 
quite large disorder. For comparison, the hopping matrix 
element in Eq.(3) is t = —D/a 2 ~ 82 meV, and numerical 
results indicate that the oscillation is almost suppressed 
when the on-site disorder to > t. The enhancement of the 
average conductance is surprising at first sight (see FIG. 
5(b)), but we attribute it to a squeezing the spatial dis- 
tribution of the edge states by disorder. Such squeezing 
can also be seen in the slight decreasing of the oscillation 
period of the conductance with large disorder strength 
as the effective radius R e ff — \/2^o7^o (see the insert 
in FIG. 5(b).). Also the disorder around the slit can in- 
crease the tunneling probability as the higher disorder 
can close the energy gap (caused by the finite size of the 
slit) near the Fermi surface [181 ]. Of course, even stronger 



disorder leads to the Anderson localization, and the edge 
states will be destroyed 0, 2^, %k > resulting in a lack of 
conductance. 

All our preceding discussion assumes strictly zero tem- 
perature. Scattering between the spin up and spin down 
edges becomes possible at nonzero temperatures due to 
inelastic scattering involving phonons, which is believed 
to cause the deviation from perfect quantization in the 
experiments of Ref. 3]. The AB oscillations discussed 
above will be observable only for disks with perimeters 
less than the phase coherence length of the edge states. 

In conclusion, we have studied theoretically a geome- 
try in which the edge states of topological insulators are 
predicted to produce coherent oscillations in the mag- 
netoconductance due to Aharonov-Bohm interference. 
This physics also produces a giant edge magnetoresis- 
tance in a weak field; FIG. 2 shows that the conductance 
may change by an order of magnitude for a tiny field 
B m 0.01T. If confirmed by experiments, this may be of 
practical interest. The giant edge magnetoresistance can 
be controlled by varying the slit width d and the size of 
the disk. We have also studied its sensitivity to disorder, 
and found that it persists for quite large disorder. The 
amplitude of AB oscillations also contains information 
about the transverse spatial extent of the edge channels. 

This work was supported by the Research Grant Coun- 
cil of Hong Kong under Grant No.: HKU 703708, and 
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